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I. INTRODUCTION 

There has been much theoretical interest in the uni- 
versal transmission fluctuations through ballistic chaotic 
systems over the past years. This activity is partially 
driven by recent experiments on electronic conductance 
in open quantum dots. Random matrix theory was shown 
to be a valuable tool to obtain analytical results on the 
distribution of transmission and reflection coefficients, as 
well as on other related quantities Q. 

Remarkably, there are very few ballistic experimental 
systems clearly showing universal transmission (or con- 
ductance) fluctuations as predicted by random-matrix 
theory. Conductance fluctuations in quantum dots 2] 
are already wanned by very small temperatures. Hence, 
theoretical clear-cut predictions of the transmission fluc- 
tuation dependence on the number of incoming and out- 
going channels 0, ^ are hardly observed. Dephasing 
effects poses further difficulties |5j , even considering that 
it can be incorporated into random matrix theory by in- 
troducing an additional phase-randomizing channel 0. 
Despite of this difficulties, quantum dots provided the 
first clear fingerprint of time-reversal symmetry breaking 
in the transmission distributions ,7]. Theory and exper- 
iment show an excellent agreement once the dephasing 
time is accounted for as a free parameter. 

An alternative to study universal transmission fluctu- 
ations is provided by microwave techniques. (There is a 
similarity to the conductante through quantum dots, that 
is proportional to the transmission - Landauer formula.) 



Transmission is directly measured in microwave experi- 
ments and cavities can be easily fabricated in any shape. 
Hence, this approach is ideally suited to verify theoretical 
predictions on transmission distributions. The first ex- 
periment of this type was performed by Doron et al. 
ft may be considered as an experimental equivalent of the 
work by Jalabert et al. on conductance fluctuations 
in essentially the same system. The first, and up to now 
only study, aiming at the channel number dependence 
and the influence of time-reversal symmetry breaking is 
our own work |lO| . For the sake of completeness we would 
like to mention that there are two further microwave ex- 
periments on non- universal aspects of transmission . 

Another quantity we shall examine in detail is the en- 
ergy derivative of the transmission, dT/dE. The mo- 
tivation stems from the study of the thermopower in 
electronic systems. There, one can show that the ther- 
mopower is proportional to the derivative of the conduc- 
tance G (or T) with respect to the Fermi energy (see, 
for instance, Ref. [Iflj for details and further references). 
Theory predicts a qualitative difference between diffusive 
and ballistic systems. Whereas for a disordered wire the 
distribution of dT/dE is expected to be a Lorentzian, for 
chaotic quantum dot systems one expects a distribution 
with a cusp at ii^ = 0. This question has been addressed 
by a number of theoretical works |0, HO, 0| . 

The comparison between random-matrix-like fluctua- 
tions and microwave experiments has limitations. It is 
not trivial to break time-reversal symmetry in microwave 
systems. On the theoretical side, on the other hand, 
analytical results are usually available for systems with 
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broken time-reversal symmetry only, whereas for systems 
with time-reversal symmetry there are formidable techni- 
cal problems. One way to break time-reversal symmetry 
in microwave systems is to introduce ferrites into the res- 
onator 0, 01 . In an externally applied magnetic field 
the electrons in the material perform a Larmor precession 
thus introducing a chirality into the system, the precon- 
dition for breaking time-reversal symmetry. It will be- 
come clear in what follows that this effect is unavoidably 
accompanied by strong absorption. 

Thus, in microwave experiments there is either no 
time-reversal symmetry breaking, or strong absorption, 
or both. Although meanwhile there is a number of works 
treating absorption 0, 0| , a better theoretical descrip- 
tion of absorption is still needed. 

Last, but not least, the coupling between the cavity 
and the waveguides is usually not perfect (or ideal) in 
the experiments. Non ideal contacts mean that part of 
the incoming flux is promptly refiected at the entrance 
of the cavity and, hence, it is not resonant. (The same 
holds for quantum dots and leads.) While most theories 
assume ideal couplin g, i t is not difficult to account for 
non-ideal coupling [llll6|. The problem, however, is that 
to quantitatively determine the quality of the contacts, 
one needs to assess the phases of the S matrix. This, in 
general, is not possible j36j. We discuss this issue in our 
analysis. 

This paper is organized as follows. In Sec. |nlwe de- 
scribe the experimental set-up and discuss how the addi- 
tion of ferrite cylinders to the microwave cavities breaks 
time-reversal symmetry. The phase-breaking features of 
the ferrite and its absorption characteristics are discussed 
in detail in App. El In Sec. IIIII we present the key el- 
ements of the statistical theory for transmission fluctu- 
ations in ballistic systems. Section Hvl is devoted to the 
statistical analysis of our experimental data. We vastly 
expand an analysis of transmission fiuctuations through 
asymmetric cavities previously presented Here we 

analyze new data on systems with reflection symmetry, 
where characteristic differences to systems without sym- 
metry are expected 0, ^3 ■ We also discuss the distri- 
bution of the derivative of the transmission with respect 
to energy, dT/dE. Our conclusions and an outlook of the 
open problems are presented in Sec. Ivl 



II. THE EXPERIMENT 

Two different cavities were used in the experiment: an 
asymmetric and a symmetric one. Refiection symmetry 
is limited by the workshop precision. Figure ^ displays 
their shapes. The height of cavities is h — 7.8mm, i.e. 
both are quasi-two-dimensional for frequencies v below 
'^max = c/2h ~ 19.2 GHz. Two commercially available 
waveguides were attached both on the entrance and the 
exit side. The cut-off frequency for the first mode is at 




FIG. 1: Sketch of the microwave cavities used in the experi- 
ments, (a) The asymmetric cavity has a = 237 mm and b can 
vary from 375 to 425 mm. (b) The symmetric one has the 
same a, while b ranges from 340 to 390 mm. The arrows indi- 
cate where the ferrite cylinders are placed. The entrance and 
exit waveguides are denoted by (1,2) and (3,4) respectively. 
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FIG. 2: Typical transmission spectrum (asymmetric cavity). 



vi — c/2w ~ 9.5 GHz where w = 15.8mm is the width 
of the wave guides. Above V2 — 18.9 GHz a second mode 
becomes propagating. All measurements have been per- 
formed in the frequency regime where there is just a sin- 
gle propagating mode. The transmission coefficients were 
measured for all four possible combinations of entrance 
and exit waveguides. Figure [3 shows a typical transmis- 
sion spectrum. By varying the length b of the resonator 
100 different spectra were taken, which were superim- 
posed to improve statistics and to eliminate non-generic 
structures. A similar procedure has been already used in 
quantum dot experiments • 

We explore the ferrite reflection properties to break 
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time-reversal symmetry: We place two hollow ferrite 
cylinders, with radius r = 10 mm and thickness d = 1 mm 
inside the cavities. The cylinders magnetization is var- 
ied by applying an external magnetic field. At an in- 
duction of B — 0.475 T the ferromagnetic resonance is 
centered at about 15.5 GHz. The electrons in the ferrite 
perform a Larmor precession about the axis of the mag- 
netic field. At the Larmor frequency the ferromagnetic 
resonance is excited giving rise to a strong microwave 
absorption. This is, of coarse, unwanted. Moving to 
frequencies located at the tails of the ferromagnetic reso- 
nance, the microwaves are partially reflected and acquire 
a phase shift depending on the sign of the propagation. 
The ferrite cylinder has thus a similar effect on the pho- 
tons as an Aharonov-Bohm flux line in a corresponding 
electron system. This correspondence has been already 
explored to study persistent currents using a microwave- 
analog [2^ . 

This method to break time-reversal symmetry has an 
obvious and unavoidable limitation: We have to move 
away from the ferromagnetic resonance frequency to 
avoid strong absorption, but have to stay close enough to 
observe a significant phase- breaking effect. In the present 
experiment, the optimal frequencies occur on a quite nar- 
row interval between 13.5 and 14.0 GHz. 

Appendix ^ gives a quantitative description of the 
phase-breaking mechanism due to the ferrite cylinders. 
Specific properties of the employed ferrite, that are use- 
ful for the understanding of the experimental data, are 
also discussed. 



III. STATISTICAL THEORY 

There are two standard statistical theories that de- 
scribe universal transmission fluctuations of ballistic sys- 
tems. One is the S'-matrix information-theoretical theory 
|23|. tailormade to calculate transmission distributions. 
The other method, where the statistical S'-matrix is ob- 
tained by modeling the scattering region by a stochas- 
tic Hamiltonian '24!|, is suited to the computation of en- 
ergy and parametric transmission correlation functions. 
Both approaches were proven to be strictly equivalent 
in certain limits ^2^1 . Complementing this result, there 
is numerical evidence supporting that the equivalence is 
general [2^. Here we use both methods: Our analytical 
results are obtained from the information-theoretical ap- 
proach, whereas the numerical simulations are based on 
the stochastic Hamiltonian one. 

We model the transmission flux deficit due to absorp- 
tion by a set of non transmitting channels coupled to 
the cavity. We consider A^i and N2 propagating modes 
at the entrance and the exit wave guides, respectively. 
The resulting scattering process is described by the block 



structured S'-matrix 

(Sii S12 Si0 \ / ~ Si0 \ 

S21 S22 S24, I = I S24, I . (1) 

s^i 5^2 s'00 J \ s^i 5^2 S00 y 

Here the set of indices {1}, {2} label the iVi, N2 propa- 
gating modes at the wave guides, while the set {(/)} labels 
the absorption channels. Transmission and reflection 
measurements, necessarily taken at the wave guides, ac- 
cess directly only the S matrix elements. 

Of particular experimental interest is the total trans- 
mission coefficient, namely, 

r-^T,b with Tab=\Sab\''. (2) 

The absorption at each channel can be quantified 
by = 1 — KS'^^)^, where (• • •) indicates an ensemble 
average (described below). We take the limits — > 00 
and 0, while keeping N^T^ = 7 constant. In this 

way we mimic the absorption processes occurring over 
the entire cavity surface, expressing their strength by a 
single parameter 7 [2^ . This modeling is equivalent to 
adding an imaginary part to the energy in the S-matrix 
|28j , a standard way to account for a finite Q- value Q . 

We obtain the distributions Pfj(T) by numerical simu- 
lation. To that end, we employ the Hamiltonian approach 
to the statistical S-matrix, namely 

S{E) = l-2T:iW\E - H + iTTWW'')-^W , (3) 

where H is taken as a member of the Gaussian orthog- 
onal (unitary) ensemble for the (broken) time-reversal 
symmetric case. This S matrix parameterization is en- 
tirely equivalent to the i^-matrix formulation recently 
used by Kogan and collaborators J^. Since the H ma- 
trix is statistically invariant under orthogonal (/3 = 1) 
or unitary (/3 = 2) transformations, the statistical prop- 
erties of S depend only on the mean resonance spacing 
A, determined by H, and the traces of W'^W . Max- 
imizing the average transmission is equivalent to put 
tr(M/'l'i^) = A/y^ ;23,]. This procedure can be used, 
in principle, to study any number N of open channels. 

The simulations are straightforward: For every real- 
ization of H we invert the propagator and compute S{E) 
for energy values close to the center of the band, = 0, 
where the level density is approximately constant. The 
dimension of H is fixed as M = 100 • • • 200, depending on 
the number of channels iV. The choice of M represents 
the compromise between having a wide energy window 
for the statistics (large M) and fast computation (small 
M). For each value of 7 we obtain very good statistics 
with 10^ • • • 10^ realizations. 

We also analyze the fluctuations of the transmis- 
sion coefficient energy derivative, dT/dE. We use the 
information-theoretical approach to analytically compute 
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moments of dT /dE. For that purpose we express dS/dE 
in terms of the S'-matrix itself and a symmetrized form 
of the Wigner-Smith time-delay matrix Qe [13, namely 



dE ~ ■ 



(4) 



Thanks to the well known statistical properties of Q^; 
matrices, the computation of {{dTab/dE)^) is possible 
[2ll |. We note that Eq. Q is strictly valid only for 
= 1. Hence, 7 is an integer number. Other values 
of 7 are obtained by extrapolation. 

The full distribution of the transmission energy deriva- 
tives, Pi3{dT/dE), is obtained by numerical simulations. 
This is a simple extension of the numerical procedure 
described above. We compute dS/dE directly from 



^ = 2mW\E 
dE 



H + inww'^y^w 



(5) 



at the same time as S{E) is calculated. 

Note that the only parameters of the theory are the 
mean resonance spacing A, the number of channels N , 
and the absorption parameter 7. In what follows we an- 
alyze the cases of asymmetric and symmetric cavities. 



Asymmetric cavities 

To this point only stochasticity and orthogonal (time- 
reversal) or unitary (broken time-reversal) symmetry 
are assumed. Additional symmetries require special S'- 
matrix parameterizations. Hence, the presented formal- 
ism is readily suited for asymmetric chaotic cavities. 

Figure El shows P{T) for the = 1 and iV = 2 cases 
for various values of the absorption 7. One can nicely ob- 
serve how the distributions for zero absorption 3] evolve 
to an exponential {N = 1) or a convolution of exponen- 
tials (iV = 2) as the absorption strength 7 increases. For 
N = 1 our simulations are in excellent agreement with 
the analytical expression obtained in Refs. [loll2^. 

For strong absorption, 7 3> 1, we find strong numeri- 
cal evidence that the distribution of individual channel- 
channel transmission energy derivatives, dTat/dE, is ex- 
ponential, namely 



P0{dTab/dE) = ^exp( 



dTab 



dE 



(6) 



where A/3 depends on 7, but not the channel indices a and 
b. Furthermore, in this regime we find that the dT„,h/dE 
for different pairs of channels are uncorrelated |31| . We 
conclude that either this distribution is insensitive to dy- 
namical channel-channel correlations, or that such corre- 
lations are insignificant in our billiards. Figure01presents 
results for typical experimental values. For independent 
dTab/dE, the distribution of dT/dE for TV = 2 is easily 
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FIG. 3: Transmission distribution P{T) for asymmetric 
chaotic cavities with N = 1 and N — 2 open channels, both 
cases with {B = 0) and without (B 7^ 0) time-reversal symme- 
try. We consider different absorption parameters 7: (solid), 
0.25 (dash), 1 (dot), 2.5 (dash dot), and 5 (dash dot dot). 
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FIG. 4: Distributions of dT/dE in units of inverse A for asym- 
metric cavities. For A'^ = 1 the distributions agree with Eq. 
@ (dotted line), while for JV = 2 they follow Eq. Q (dotted 
line). 



obtained by a convolution using Eq. © and reads 
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dT 
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(7) 



It remains to relate A/3 to 7. This is done by computing 
{{dTab/dEY). The latter can be analytically calculated 
using the energy derivative of the S'-matrix, Eq. Q , and 
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reads l31 



dE 



TT 



A2 Q;2(a + 1)2 

a2 + a - 2 + 4(2 - /?) 



a 



4(2 - 13) 



(8) 



where a = TVi + iV2 + 7. Recalling that {{dTab/dEf) = 
2/A| we find as a function of 7. 

In Fig. 01 we compare the approximation Pp{dT / dE), 
where \p calculated as described above, with a direct 
numerical simulation. The agreement is rather good. 



Symmetric cavities 
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The influence of absorption on the transmission fluctu- 
ations is even more pronounced in billiards with reflection 
symmetry. In the absence of absorption the transmis- 
sion distributions for reflection symmetric cavities was 
already analytically computed. The most salient fea- 
tures are the following: When time-reversal symmetry is 
preserved, the theory predicts that the transmission dis- 
tribution P(T) for reflection-symmetric cavities remains 
invariant when T is substituted by 1 — T [s^- On the 
other hand, for broken time-reversal symmetry, P(T) co- 
incides with the one for the asymmetric case, but with T 
replaced by 1-T 18]. 

To fulfill the reflection symmetry, it is sufficient to con- 
sider the S'-niatrix with the block structure 1331 



(9) 



where Si and S2 are unitary (and symmetric for f3 — 1) 
Nt/2 X Nt/2 matrices with Nt = 2iV -I- iV^. Both Si 
and 5*2 have the structure given by Eq. l(T]l. 
The transmission coefficient now reads 



^ — T S |['S'l]afc — [<5'2]afc| = ^ CTa, 



(10) 



„b=l 



We numerically generate P/iiT) and Pp{dT / dE) using 
the Hamiltonian approach to the S-matrix, Eqs. ^ and 
Q. Now two statistically independent matrices. Si and 
52, are required. We chose the dimension of H to be AI = 
50. For each value of 7 we obtain very good statistics 
with 10^ realizations. 

FigureElcontrasts Pp{T) obtained analytically for zero 
absorption 0| with our numerical simulations for differ- 
ent values of 7. Our analysis is restricted to the N = 1 
and iV = 2, as before. We observe that with increasing 
7 the fingerprints of the reflection symmetry fade away, 
and the distributions become quite similar to those of 
asymmetric cavities. 



FIG. 5: Symmetric cavity transmission distributions P(T) 
for the one- and two-channel case. For B = we consider 
7 = 0, 0.5, 2, 2.5, and 10, corresponding to the solid, dashed, 
dotted, dashed-dotted, and dashed-dotted-dotted lines re- 
spectively. The same for the case of broken time-reversal 
symmetry, B 0, but with 7 = 0, 0.25, 1, 2.5, and 5. 



As in the asymmetric case, for the strong absorp- 
tion regime, 7 3> 1, our numerical simulations strongly 
suggest that the distribution of the energy derivative 
of individual channel-channel transmission coefficients 
Pi3{dTab/dE) is exponential. However, in distinction to 
the asymmetric case, here the exponential law depends 
on the channels: The reflection symmetry (see Fig. ^ 
makes the channels (1,4) and (2,3) indistinguishable. Ac- 
cordingly, we find that the "diagonal" coefficients T14 and 
T23, denoted by crjjj, and the "off-diagonal" ones T24 and 
T13, denoted by (7°^ have different variance. The second 
moment of the diagonal da'^f^/dE is [sj 



V dE 



A2 (a'-2)a'^(a'-f 1)2 



a'(a'-l)(7/?-6) (a' V a' + 2)(2 - /3) 



whereas the off-diagonal is 



1 



2{a'^ + a' + 2) 



(11) 



dE 



A2 (a'-2)a'2(a'-f l)2(a' + 3) 

(2-/3)^-f4(/?-l)l. (12) 
a' -I- 1 



Here, a' = /3(iV-F7/2). 

For 7 3> 1, based on the numerical simulations, we 
assume that P[}{daab/dE) is exponential and that for 
different pair of channels a and b the duab/dE are un- 
corrected. We then equate ^| = 2/{{da^f^/dE)^) and 
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FIG. 6: Transmission energy derivative distributions for the 
symmetric chaotic cavities. The points represent the results 
of the simulations for 7 = 18 (22) for B = (B / 0) for 
= 1; 7 = 14 (18) for B = (B / 0) for N = 2. The dotted 
lines give the approximations © and 1131 . For = 1 we 
present the diagonal case. 



i^l = 2/{{d<Tl^/dEf) to write 



Pp{dT/dE) 
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1 1 \ ( 

1 exp --^ 



dT 



dE 
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dT 


ai 


"2/ 




a2 




dE 



X exp 





dT 






dE 


)] 



(13) 



where ai = + a2 = M/3 ~ ^p/'^- Figure El com- 

pares the approximation Pp{dT / dE) with our numerical 
simulations. We chose parameters realistic to out experi- 
ment. The agreement is quite good. Deviations between 
the approximation and the numerical simulations 
are of order I/7. 



IV. STATISTICAL ANALYSIS OF THE 
EXPERIMENTAL RESULTS 



The statistical analysis of our experiment is based on 
two central hypothesis. First, as standard, we assume 
that the transmission fluctuations of a chaotic system are 
the same as those predicted by the random matrix theory 
[2^. Second, we employ an ergodic hypothesis to justify 
that ensemble averages are equivalent to running aver- 
ages, that is, averages over the energy (frequency) and/qr 
shape parameters. This requires RMT to be ergodic '3^, 
which was recently shown |39j. With few exceptions (see, 
for instance, Ref. [s^]) this point is unnoticed. 

The experimental transmission coefficients were ob- 
tained by superimposing 100 different spectra measured 



FIG. 7: Mean Transmission (T) for the TV = 2 case for B = 
(sohd line) and B = 0.470 T (dotted line). The Larmor 
resonance frequency is uir = 2-k x 14.86 GHz) 



for billiard lengths b (see Fig. In the studied fre- 
quency regime there is only a single propagating mode 
in each of the waveguides. Hence, to every waveguide 
we associate a single scattering channel. For the = 1 
case all measurements for the different combinations of 
entrance and exit waveguides were superimposed. The 
transmission for the N = 2 case was obtained by com- 
bining the results from all = 1 measurements, namely, 
T = Ti3 + Ti4 + T23 + T24,. 

Figure [7| shows the mean transmission {N — 2 case) 
with and without applied external magnetic field. When 
related to experimental quantities, (• • •) indicate running 
averages. Using the Weyl formula, we associate the fre- 
quency V (actually, i^/A) with the energy E introduced 
in the preceding section. The strong absorption due to 
the Larmor resonance is clearly seen. In App. ^ we 
discuss why is the phase-breaking effect expected to be 
best observed in the tails of the Larmor resonance. Fig- 
ure IHl illustrates this very nicely. It shows the scaled 
transmission distribution P(T/{T)) for the asymmetric 
billiard in three different frequency windows both with 
and without applied external magnetic field. It is only 
in the frequency interval form 13.55 to 13.85 GHz that 
P{T/{T)) changes with magnetic field. We stress that 
this is different from just an absorption effect: In the 
frequency window around 14.45 GHz, where the absorp- 
tion is strongest, the normalized distributions with and 
without magnetic field are basically the same (the only 
difference is in the mean transmission). We identify the 
change in P{T/{T)) with the expected phase-breaking 
effect and assume that the applied magnetic field is suffi- 
cient for the ferrite cylinders to fully break time-reversal 
symmetry. Similar observations were made for the sym- 
metric billiard. 

Before we present our statistical analysis, it remains 
to discuss how ideal the cavity-waveguides coupling is. 
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FIG. 9: Transmission distributions for the asymmetric cavity. 
The histograms correspond to data taken within the indicated 
frequency window. The dotted lines stand for the random 
matrix simulations, with 7 as a fitting parameter. 



FIG. 8; Transmission distribution for the A'' = 1 (left) and 
N = 2 (right) channel cases for three different frequency win- 
dows of width Si/ — 0.3 GHz centered at vo (indicated in the 
figure). The histograms correspond to B = (solid line) and 
B = 0.470 T (dotted line). 



To determine the antenna coupling we measured the 
transmission through two waveguides facing each other 
directly. In the whole applied frequency range the to- 
tal transmission was unit, with an experimental uncer- 
tainty below 5%, showing that the antenna coupling is 
perfect. There are, however reflections of about 10% in 
amplitude from the open ends of the waveguides, where 
they are attached to the billiard. Small deviations from 
ideal coupling are also consistent, for the frequencies we 
work, with Ref. Is^ . Since the absorption is strong in 
the present experiment, and an imperfect coupling can 
be compensated for to a large extent by a rescaled ab- 
sorption constant, we decide not to explicitly account for 
coupling corrections. In summary, throughout the forth- 
coming analysis we assume perfect coupling between the 
cavity and the waveguides. 

For the sake of clarity, we present the statistical anal- 
ysis of the asymmetric and the symmetric cavities sepa- 
rately. 



Asymmetric cavity distributions 

Figure 1^ compares the experimental transmission dis- 
tributions in the "phase-breaking" frequency window 
with the statistical theory. The absorption parameter 
7, see Sec. IIIII was adjusted to give the best fit of the 
theoretical (T) to the experiment. The agreement is ex- 
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FIG. 10: Distribution of the energy derivative of the transmis- 
sion for the asymmetric cavity. The dotted lines correspond 
to the theoretical distributions. 



cellent, except for = 2 with B ^ 0. 

We work with a single asymmetric cavity, but use dif- 
ferent 7 values for = 1 and N = 2. The reason is 
simple: For N = 2 we consider the contributions from all 
antennas to the transmission, whereas for iV = 1 two an- 
tennas act as additional absorption channels. This gives 



{N=l) 



7 



{N=2) 



rise to a simple relation, namely, 7 

In order to compare the experimental transmission 
energy derivative distributions with the universal ran- 
dom matrix results we have to rescale the experimental 
data by the mean resonance spacing, namely, dT/dE 
AdT/dE. We use the Weyl formula to estimate A. Fig- 
ure 1101 shows a comparison between theoretical and ex- 
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^= 1 13.55 - 13.85 GHz 





N=2 13.55- 13.85 GHz 




FIG. 11: Normalized joint distribution F{T,dT/dE) = 
P{T, dT/dE)/[P{T)P{dT/dE)] for the asymmetric cavity for 
N = 1, B = Q. Similar result holds to B 7^ 0. 



FIG. 13: Transmission distributions P{T) for the symmetric 
billiard. Histograms stand for the data taken at the indicated 
frequency interval, whereas the dotted lines correspond to the 
simulations. The absorption 7 is a fitting parameter. 




FIG. 12: Same as in Fig. [TT] but with dE/dT replaced by 
{dT/dE)resc = {dT/dE)/^T{l-T). 

perimental results for P{AdT/dE). Note that we take 
the same 7 as for P{T). The signatures of the chan- 
nel number, and the influence of time-reversal symmetry 
breaking are clearly seen. We checked that the increase in 
absorption when switching on the magnetic field, without 
switching to the unitary ensemble as well, is not sufficient 
to reproduce the data. Inaccuracies in the assessment of 
A provide a possible explanation for the slight disagree- 
ment between theory and experiment. The Weyl formula 
does not account for the standing waves in the ferrite 
cylinders and, thus, overestimates A. This is consistent 
with Fig. Cni 

The joint distribution of T and dT/dE was studied 
in Ref. for TV = 1 and 7 = 0. Remarkably, it 
was found that albeit T and dT/dE are correlated, the 
rescaled quantity {dT/dE)rcsc = (dT/dE)/ ,/T{l - T) 
and T are not. We checked if this finding holds in 



our experiment, despite of absorption. Figure ITTI shows 
the "normalized" joint probability F{T,MT/dE) = 
P{T,AdT/dE)/[P{T)P{AdT/dE)] in a three dimen- 
sional representation for = 1 and B ^ 0. A clear 
correlation is observed. To contrast. Fig. El shows 
F[T, A{dT / dE)rcsc]- Here the distribution becomes flat. 
Unfortunately we do not have enough statistics to make 
a reliable determination of the distribution. A similar 
result, not shown here, holds for the B case. 



Symmetric cavity distributions 

We switch now to the statistical analysis of the sym- 
metric cavity transmission fluctuations. 

Figure E| shows the experimental shows P{T) for 
transmissions within 13.55 < i' < 13.85 GHz, where the 
phase-breaking effect is expected to be strongest. As be- 
fore, the absorption parameter 7 is the best fit of the the- 
ory to the experiment. Here, for all studied cases a nearly 
perfect agreement is found. Now 7(^=1) = ry(N=2) ^ 
This is due the reflection symmetry. 

Figure 1141 shows the experimental distributions 
P{AdT/dE) for the symmetric case. The signatures of 
the channel number and the influence of breaking time- 
reversal symmetry, are clearly seen. For all cases of the 
symmetric billiard the theoretical curves are plotted as 
well. We observe that the experimental distributions ver- 
ify the overall trends of the theoretical predictions. In 
particular, the characteristic cusp at £' = is nicely re- 
produced for = 1 . Similar to the asymmetric case, the 
agreement between experiment and theory is not as good 
as for the transmission distribution. 

As in the case of asymmetric cavities, theoretical 
calculations show that although T and dT/dE 
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13.55 - 13.85 GHz 




-0.1 0.0 0.1 

MT/dE 



0.2 



0.15 0.00 0.15 0.30 

AdT/dE 




FIG. 14: Distribution of the energy derivative of the transmis- 
sion for the symmetric cavity. The dotted hnes corresponds 
to the theoretical distributions obtained from random matrix 
theory. 



FIG. 16: Same as in Fig. [T^ but with dE/dT replaced by 
{dT/dE)resc = {dT/dE)/^T{l-T). 



V. CONCLUSIONS 




FIG. 15: Normalized joint distribution F{T,dT/dE) = 
P{T,dT/dE)/[P{T)P{dT/dE)] for the symmetric cavity for 
iV = 1, B = 0. Similar resuh holds to B / 0. 



are correlated, the rescaled quantity {dT / dE)^csc = 
{dT/dE)/ .JT{1 - T) is independent of T. Here also 
the analytical results were obtained for the iV ^ 1 
case. Fig. 1151 shows the normalized joint probabil- 
ity F(T,dT/dE) = P{T,dT/dE)/[P{T)P{dT/dE)] in a 
three dimensional representation for N = 1, B = case. 
A clear correlation is manifest. For comparison, Fig. 1161 
shows the corresponding quantity for f[T, {dT/dE)^csc]- 
Now the correlation has vanished, in accordance with the- 
ory. Similar result, not shown here, holds for the B ^ 
case. 



This work shows that microwaves are ideally suited to 
experimentally verify the theory of universal transmis- 
sion fluctuations through chaotic cavities. The results 
presented in the present paper would have been hardly 
accessible by any other method. 

We observe a nice overall agreement between our ex- 
perimental data and the random matrix results. How- 
ever, the comparison between theory and microwave ex- 
periment is limited by the following issues. 

In experiments, the coupling between waveguides and 
the cavity is usually not ideal, whereas in most theoret- 
ical works ideal coupling is assumed. In the frequency 
range we work (36| supports our working hypothesis of 
nearly perfect coupling. In general, however, it turns out 
that without measuring the S matrix (with phases) it 
is hard to disentangle direct reflection at the cavity en- 
trance (imperfect coupling) from absorption. From the 
experimental side, it would be desirable to have a better 
handle on absorption. 

Microwave systems are usually time-reversal invari- 
ant, and as we have seen it is not trivial to break this 
symmetry. At the same time we increase the magnetic 
field, turning on the phase-breaking mechanism, absorp- 
tion also increases. Unfortunately, both effects are in- 
extricable. This is why it is beyond our present experi- 
mental capability to quantitatively investigate the trans- 
mission fiuctuations along the crossover regime between 
preserved and broken time-reversal invariance. Actu- 
ally, to compare theory with experimental results we as- 
sume that the transmission data at B = 0.470 mT and 
13.55 < ly < 13.85 GHz are far beyond the crossover 
regime. 

We hope that the present work will trigger additional 
theoretical effort in the mentioned directions. 
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APPENDIX A: PHASE-BREAKING 
PROPERTIES OF THE FERRITE 

This appendix is devoted to the discussion of the ferro- 
magnetic resonance and the phase-breaking mechanism. 
For that purpose we first quickly present some elements 
of the well-established theory of microwave ferrites, see 
for instance, Ref. [s^ 

For the sake of simplicity, we first restrict ourselves to 
the situation of an incoming plane wave reflected by the 
surface of an semi-infinite ferrite medium. We assume 
that incoming, reflected, and refracted waves propagate 
in the xy plane and are polarized along the z direction, 
and that there is an externally applied static magnetic 
field in the z direction, as shown in Fig. El We ask what 
is the phase acquired due to the reflection on the ferrite. 

To answer this question we need to solve Maxwell's 
equations. For this geometry and single-frequency elec- 
tromagnetic fields, like our microwaves, this is a simple 
task. The ferrite properties come into play by the consti- 
tutive relations D = eoeE and B — /lo/iH, more specif- 
ically through the permeability /x, that is a tensor with 
the form 



1+Xr 

^Xl 



with 



Xr = 



+ Xr 



Xo 



uj ^u + iX. (A2) 



Here ui^ = —jHq and ujm — jMg are the precession 
angular frequencies about the external field Hq and the 
equilibrium magnetization Mq, respectively, /io is the 
static susceptibility. More details can be found, for in- 
stance, in Chapter 2.2.3 of Ref. i34i] . 

We solve the proposed problem using for the electric 
field the ansatz E(r) = E{r)e2,, where 



X <0, 
x> 0, 



(A3) 



with k/ = ko{— cos a, sin a , 0) , k/f = fco(cos a, sin a , 0) , 
kr = A;(— cos/3, sin/3,0), see Fig.lTTI 

The derivation of the amplitudes Ej.Er and Et is sim- 
ilar to that of Fresnel's formula (see, for instance, 1351) . 
Since an explicit calculation for ferrites is given in |22l |. 




* X 



FIG. 17; Plane wave reflected by the surface of a ferrite slab. 



only the results shall be given. Using the continuity of 
E[|, Dj^, and H[| on the boundary, one writes 

Et = Ei + Ell and fc sin /? = fco sin a (A4) 

which is just Snell's law. For the relative amplitude of 
the reflected part we obtain 



Er 
Ei 

where 



(n^/e) cos a + iS sin a — \/ — 



sin^ OL 



[rP- 1 e) cos a — i6 sin a -I- 



sin^ a. 



n 



(Al) and 



(5 = 



1+Xr LUl{lOl+L^m) 



(A5) 



(A6) 



(A7) 



Note that there is a term depending on the sign of a, 
i.e. on the direction of the incident wave. This term is 
responsible for the phase-breaking effect. 

The above formulas have to be modified when dealing 
with a ferrite of finite width. For a slab of thickness I 
and a — Q we have 



Ej" 



4-^ 



(1 ^^2gife(l-n)i — (1 — i.)2gife(l-|-ra)i 



(A8) 



and 

Er 
Ei 



-2i sin knl- 



1 



' (1 _|_ e.)2gife(l~ri)i — (1 — .l)2gife(l+ri)i " 

(A9) 

In contrast to Eq. (|A5|I . Et is no longer the amplitude of 
the transmitted wave propagating inside the ferrite. Here 
Et is the amplitude of the wave that crossed the ferrite 
slab an emerged at the other side. The explicit formula 
for a 7^ is lengthy and is not be presented here. 
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12 14 16 

Frequency p in GHz 



FIG. 18: Reflection, transmission, and phase shift for a ferrite 
slab (Mo =130 mT, e =15, A = 0.1 GHz) of thickness I =1 
mm at Bo = Mo-f^o =470 mT for different incidence angles a. 



The phase-breaking becomes clearly manifest by writ- 
ing Eq. HA9|I as 



Er 
Ei 



Er 
Ei 



i0rotl(a). 



(AlO) 



where 0rcfl(Q!) is the phase acquired due to reflection. 
Figure ITHl shows modulus of transmission \Et/Ei\ and 
reflection \Er/Ei\ as well as the phase shift for differ- 
ent incidence angles and I — 1 mm, the thickness of 
our ferrite cylinders. The curves are calculated using 
the ferrite parameters (see caption of Fig. I18|l given by 
the supplier. We find a resonance angular frequency of 
ljJr = ^ujl{uol + ljJm) ~ 271 X 14.86 GHz. This reso- 
nance corresponds to the dominant structure observed in 
Fig.^1 The additional substructures are due to standing 
waves inside the ferrite. 

To illustrate the phase-breaking effect of the ferrite, in 
Fig. El we show the phase diflFerence Acj) = (j)^efi{ct) — 
(/)refl(— a) between the incoming and the time-reversed 
wave. We see that the effect is maximal at the resonance 
frequency, and vanishes as one moves off-resonance. Un- 
fortunately, the absorption is maximal at the resonance 
too. This are the quantitative observations in support of 
the discussion presented in Sec. m 
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FIG. 19: Difference Acp = (;/>rcfl(a) — <^rofl(— q) of phase shifts 
observed between an incoming wave and its time-reversed 
equivalent. 



(a) 




12 14 16 

Frequency i/ In GHz 



(b) 




12 14 16 
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FIG. 20: Experimental reflection for a ferrite slab of thick- 
ness I =lmm (a) and 2 mm (b). The dashed lines have been 
calculated by superimposing the results for two different in- 
ternal magnetizations A'lo =110 mT an 190 mT. The broad 
minimum observed for / =2 mm close to 13 GHz is due to 
a standing wave within the ferrite. For / =1 mm the corre- 
sponding minimum is at 15 GHz and superimposes the ferro- 
magnetic resonance. 
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Finally, to experimentally check the properties of the 
ferrites, we place a small sheet of the material between 
two waveguide facing each other. Two different thick- 
nesses I = 1 mm and 2 mm were used. Figure |^ shows 
the measured reflection \Efi/Ej\ as a function of i'. 
The small oscillations superimposing the dominant res- 
onance structures correspond to standing waves within 
the waveguide and are an artifact of the experiment. 
Comparing the experimental results with the calculation 
shown in Fig. 1181 we notice that the assumption of a 
single homogenous internal magnetization is not in accor- 
dance with the measurement. The dashed line is obtained 
by superimposing the theoretical results for two different 
values of the magnetization. The overall behavior of the 
resonance structures becomes then in qualitatively agree- 
ment with the data. 
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